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Analytical solutions are presented for a class of generalized r— modes of rigidly rotating uniform 
density stars — the Maclaurin spheroids — with arbitrary values of the angular velocity. Our analysis 
is based on the work of Bryan M ; however, we derive the solutions using slightly different coordinates 
that give purely real representations of the r— modes. The class of generalized r— modes is much 
larger than the previously studied 'classical' r— modes. In particular, for each I and m we find l — m 
(or I — 1 for the m = case) distinct r— modes. Many of these previously unstudied r— modes (about 
30% of those examined) are subject to a secular instability driven by gravitational radiation. The 
eigenfunctions of the 'classical' r— modes, the I — m + 1 case here, are found to have particularly 
simple analytical representations. These r— modes provide an interesting mathematical example of 
solutions to a hyperbolic eigenvalue problem. 

PACS Numbers: 97.10.Sj, 97.10.Kc, 04.30.Db 
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I. INTRODUCTION 

During the past year the r— modes of rotating neu- 
tron stars have been found to play in interesting and 
important role in relativistic astrophysics. Andersson |l]] 
and Friedman and Morsink || showed that these modes 
would be driven unstable by gravitational radiation re- 
action in the absence of internal fluid dissipation. Lind- 
blom, Owen, and Morsink || have subsequently shown 
that this instability will in fact play an important role in 
the evolution of hot young neutron stars. The gravita- 
tional radiation reaction force in these modes was shown 
to be sufficiently strong to overcome the internal fluid dis- 
sipation present in neutron stars hotter than about 10 9 K. 
Hot young rapidly rotating neutron stars are expected 
therefore to radiate away most (i.e. up to about 90%) of 
their angular momentum via gravitational radiation in a 
period of about one year. Owen, et al. 0] have shown 
that the gravitational radiation emitted during this spin- 
down process is expected to be one of the more promising 
potential sources for the ground based laser interferom- 
eter gravitational wave detectors (e.g. LIGO, VIRGO, 
etc.) now under construction. 

To date the various analyses of the r— modes and their 
instability to gravitational radiation reaction have all 
been based on small angular velocity approximations. 
This instability is of primary importance in astrophysics 
for rapidly rotating stars. The purpose of this paper is to 
provide the first look at the properties of these important 
modes in stars of large angular velocity. We do this by 
solving the stellar pulsation equations for the r— modes 
of the rapidly rotating uniform density stellar models 
which are known as the Maclaurin spheroids. The pul- 
sations of these models were studied over a century ago 
by Bryan |^] , who showed how analytical expressions for 



all of the modes of these stars could be found. We follow 
the general strategy developed by Bryan to derive analyt- 
ical expressions for the r— modes of these stars. We use 
somewhat different coordinates than Bryan, however, in 
order to obtain real representations of the r— modes (of 
primary interest to us here) using purely real coordinates. 

We generalize the traditional definition of r— mode to 
include any mode whose frequency vanishes linearly with 
the angular velocity of the star. Such modes have as their 
principal restoring force the Coriolos force, and hence it 
is appropriate to call them rotation modes or generalized 
r— modes ||. We find a very large number of generalized 
r— modes in the Maclaurin spheroids In particular 
for each pair of integers I and m (with I > m > 0) 
we find I — m (or I — 1 for the m = case) distinct 
r— modes. The 'classical' r— modes as studied for exam- 
ple by Papaloizou and Pringle M correspond here to the 
case I = m + 1 ||] . We find that many of the previously 
unstudied r— modes (about 30% of those examined) are 
also subject to the gravitational radiation driven instabil- 
ity in stars without internal fluid dissipation. These new 
r— modes couple to higher order gravitational multipoles 
and consequently are expected to be of less astrophysi- 
cal importance than the I = m + 1 = 3 mode that is of 
primary importance in the instabilty discussed by Lind- 
blom, Owen, and Morsink ||. 

In Section [n] we review a few important facts about 
the equilibrium structures of the Maclaurin spheroids. 
In Sec. |H we present the equations for the modes of 
rapidly rotating stars using the two-potential formal- 
ism p0| . This formalism determines all of the proper- 
ties of the modes of rotating stars from a pair of scalar 
potentials: a hydrodynamic potential SU and the grav- 
itational potential The equations satisfied by these 
potentials are (coupled) second-order partial differential 
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equations with suitable boundary conditions at the sur- 
face of the star and at infinity. These equations be- 
come extremely simple in the case of uniformly rotating 

we introduce coordi- 



IV 



uniform-density stars. In Sec. 
nates which allow the equations for the two potentials 
to be solved analytically. And following in the footsteps 
of Bryan ||, we present the general solutions to these 
equations for the generalized r— modes of the Maclaurin 
spheroids. In Sec. [v] we give expressions for the various 
boundary conditions in the coordinates adapted to this 
problem. In Sec. [v| we deduce the eigenvalue equation 
that determines when the boundary conditions may be 
satisfied. In Sec. [Vl| we also present explicit solutions to 
this eigenvalue equation for a large number of general- 
ized r— modes. In the limit of small angular velocity we 
tabulate a complete set of solutions for all generalized 
r— modes with I < 6. We also present graphically the 
angular velocity dependence of each of these r— modes 
which is unstable to the gravitational radiation instabil- 
ity. In Sec. VII we analyze the analytical expressions for 



the eigenfunctions of these modes. We show that in the 
case of the 'classical' r— modes, the I = m + 1 case here, 
the eigenfunctions 5U and d<& have particularly simple 
forms. In particular, each of these eigenfunctions is sim- 
ply r m+1 Y m+ i m (9, ip) (where r, 9 and ip are the standard 
spherical coordinates) multiplied by some angular veloc- 
ity dependent normalization. We also show that these 
'classical' r— mode eigenfunctions have the unexpected 
property that Ap, the Lagrangian pressure perturbation, 



vanishes identically throughout the star. In Sec. VIII 



discuss some of the interesting implications of this anal- 
ysis. And finally, in the Appendix we explore in some 
detail the properties of the rather unusual bi-spheroidal 
coordinate system needed in Sec. IV to solve the pulsa- 
tion equations for the hydrodynamic potential 5U . 



II. THE MACLAURIN SPHEROIDS 

The uniformly rotating uniform-density equilibrium 
stellar models are called Maclaurin spheroids. The struc- 
tures of these stars are determined by solving the time 
independent Euler equation 



= V a 



i(^+j/ 2 )0 2 + ^-$ 



(2.1) 



In this equation p is the pressure, p is the density, f2 is the 
angular velocity, and $ is the gravitational potential of 
the equilibrium star. Using the expression for the grav- 
itational potential of a uniform-density spheroid pill, it 



is straightforward to show that the solution to Eq. (2.1) 
for the pressure is 



p = 2nGp 2 C(l + C 2 )(l - Co cot" 1 Co) 



x a 



x + y z 

T+cT~cI 



(2.2) 



where a is the focal radius of the spheroid, G is Newton's 
constant, and Co is related to the eccentricity e of the 
spheroid by e 2 = 1/(1 + Co)- Similarly, it follows that the 
angular velocity of the star is related to the shape of the 
spheroid by 



fl 2 = 2nGp( [(1 + 3C 2 ) cot" 1 Co - 3Co] 



(2.3) 



We note that small angular velocities, O, correspond to 
small eccentricities e and large Co- 

The surfaces of these stellar models are the surfaces on 
which the pressure vanishes: 



y 



C 2 + i C 2 



= a 2 . 



(2.4) 



This surface is an oblate spheroid. Let us denote the 
equatorial and polar radii of t his s pheroid as R e and R p 
respectively. We see from Eq. (|2.4[) that 



R,, 



« 2 (Co 2 + l), 



R p 



a 2 C 



(2.5) 
(2.6) 



If we consider a sequence of uniformly rotating spheroids 
having the same total mass, then the volume of each 
of the spheroids in the sequence is the same (since the 
density is constant). Let R denote the average radius of 
the spheroid: R 3 = R\R p . It follows that R is constant 
along this sequence since the volume of a spheroid is V — 
AttR 3, /3. Thus the angular velocity dependence of the 
focal radius a is determined by 



R 4 



(2.7) 



since Co is rela ted to the angular velocity of the spheroid 
£1 by Eq. ( |2.3| ). This expression determines then the an- 
gular velocity dependencies of the equatorial and polar 
radii of the spheroid: 



R e — R 



Rp — R 



Co 2 + l 
Co 2 



C 2 + i 



1/6 



1/3 



(2.8) 



(2.9) 



In the work that follows we will need the quantity 
rt a V a p, where n a is the outward directed unit normal 
to the surface of the star, in order to evaluate certain 
boundary conditions associated with the stellar pulsa- 
tions. Since V a p is also normal to the surface of the 
star, we may use the expression (n a \7 a p) 2 = S7 a p\7 a p 
and Eq. dO) to obtain 



(2.10) 



n a V aP = -4^Gp 2 C 2 (l+Co 2 )(l " Co cot- 1 Co) 



(1 



H/ 2 z 2 
C 2 ) 2 Co 4 



1/2 



where ( x, y, z) are to be confined to the surface defined 
by Eq. (Oh. 
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III. THE PULSATION EQUATIONS 



K 2 V a X7 a SU - Az a z b V a V b SU = 0, 



(3.6) 



The modes of any uniformly rotating barotropic stellar 
model are determined completely in terms of two scalar 
potentials SU and (5$ Q . The potential S<& is the New- 
tonian gravitational potential, while SU is a potential 
that primarily describes the hydrodynamic perturbations 
of the star: 



SU 



Sp 

p 



5<f>, 



(3.1) 



where Sp is the Eulerian pressure perturbation, and p is 
the unperturbed density of the equilibrium stellar model. 
We assume here that the time dependence of the mode is 
e tuJt and that its azimuthal angular dependence is e mv , 
where to is the frequency of the mode and m is an integer. 
The velocity perturbation 5v a is determined by solving 
Euler's equation. This reduces in this case to an algebraic 
relationship between Sv a and the potential SU : 



Sv a = iQ ab V b SU. 



(3.2) 



The tensor Q ab depends on the frequency of the mode, 
and the angular velocity of the equilibrium star f2: 



Q 



ab 



i 



(uj + m n) 2 - An 2 
(lu + m n)s ab - 



An 2 



-z a z b 



2iV a v b 



(3.3) 



In Equation (3^3) the unit vector z a points along the ro- 
tation axis of the equilibrium star, S ab is the Euclidean 
metric tensor (the identity matrix in Cartesian coordi- 
nates), and v a is the velocity of the equilibrium stellar 
model. The potentials SU and S$ are determined then 
by solving the perturbed mass conservation and gravita- 
tional potential equations. In this case, these reduce to 
the following system of partial differential equations [ |l~C| ] 

V a (pQ ab VbSU) = -(u + mn)p^f-(SU + 5$), (3.4) 

dp 



V a V a £$ = -4nGp^f-(SU + 6$). 

dp 



(3.5) 



In addition these potentials are subject to the appropri- 
ate boundary conditions at the surface of the star for SU 
and at infinity for (5$. 



The stellar pulsation Eqs. (3^4) and (3__5) simplify con- 
siderably for the case of uniformly rotating uniform- 
density stellar models. In this case dp/ dp = pS(p) [where 
8(p) i s th e Dira c de lta function]. Thus the right sides of 
Eqs. (3.4) and ( |3.5| ) vanish except on the surface of the 
star. Further, the density p that appears on the left side 
of Eq. (3.4) may be factored out. The resulting equations 
then in the uniform-density case are simply 



V Q V a 6$ = -4nGp 2 S(p) (SU + <5$) , (3.7) 

where k is related to the frequency of the mode by 

«f2 = to + mfi. (3-8) 

These equations are equivalent to those used by Bryan || 
in his analysis of the oscillations of the Maclaurin 
spheroids. 

Next we wish to consider the boundary conditions to 
which the functions SU and <5$ are subject. Let £ denote 
a function which vanishes on the surface of the star, and 
which has been normalized so that its gradient, n a — 
V a S, is the outward directed unit normal vector there, 
n a n a = 1. First, the function SU must be constrained 
at the surface of the star, £ = 0, in such a way that the 
Lagrangian perturbation in the pressure vanishes there, 
Ap = 0. This condition can be written in terms of the 
variables used here by noting that 



Ap 



, Sv a 



(3.9) 



Then using Eqs. ( p.l| ) and (3.2) the boundary condition 
can be written in terms of SU and <5$ as 



= 



pnVL(5U + 5$) + Q ab V aP V b SU 



(3.10) 



S|0 



The perturbed gravitational potential 5<& must vanish 
at infinity, lim r ^ 00 5$ = 0. In addition (5$ must as a 



consequence of Eq. (3.7) have a finite discontinuity in its 
first derivative at the surface of the star. In particular 
the derivatives must satisfy 



n a V a S<f> 






EJ.0 



7i a V a <S$ 



AirGp 2 {8U + <5$) 



n a V a p 



STO 
(3.11) 



The problem of finding the modes of uniform-density 
stars is r educe d th erefore to finding the solutions to 
Eqs. (3.6) and ( ]3.7| ) subject to the boundary conditions 
in Eqs. Q3.1C|) and ( pl| ). 



IV. SOLVING FOR THE POTENTIALS 

In this section we find the general so luti ons fo r th e 
potentials SU and 5<& that satisfy Eqs. (3J3) and ([3.7]). 
This analysis basically follows that of Bryan H except 
for some changes to modernize notation, and a change 
of coordinates to express in a purely real manner the 
solutions of interest to us here. Our primary concern 
here is to find expressions for the generalized r— modes 
of the Maclaurin spheroids. 
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We first introduce a system of spheroidal coordinates 
that are useful in solving for the gravitational potential 
<5<f>. Thus we introduce the coordinates (/z, £, <p>) that are 
related to the usual Cartesian coordinates (x, y, z) by the 
transformation: 



x = ay/ (C 2 + 1)(1 - n 2 ) cos(p, 



y = aV(C 2 + l)(l-M 2 ) 



z = aQu, 



sin if, 



(4.1) 



(4.2) 



(4.3) 



where a is defined in Eq. (2^) above. These are the stan- 
dard oblate spheroidal coordinates. It is straightforward 
to show that 



y 



i + C 2 C 2 



a 2 . 



(4.4) 



Thus the surfaces of constant C are oblate spheroids, with 
C = Co corresponding to the surface of the star. The 
coordinate C has the range < C < Co within the star, 
and C > Co outside. The surface C = corresponds to a 
disk of radius a within the equatorial plane of the star. 
The nature of the surfaces of constant /x can similarly be 
explored by noting that 



Q" l (iC) are bounded in the limit C — > 0. However, we 
must insure that the solution is smooth across the disk 
C = 0. The functions Qf l (iQ are non-zero at iC = 
(see Bateman |§ Eqs. 3.4.9, 3.4.20 and 3.4.21). For the 
case of l + m odd the function Q^^QP^in) is therefore 
discontinuous at the disk C = 0, and consequently it does 
not satisfy Laplace's equation there. Similarly for I + m 
even the function Q™(iC)P™(/i) is continuous but not 
differentiable at C = 0, and so it does not satisfy Laplace's 
equation there either. Thus, in the interior of the star, 
C < Co, the solution to Eq. (3.7) for given I and m is 



<5$ 



~pr«o) 



(4.8) 



The potentials in Eqs. (4.7) and (4.8) have been normal- 
ized so that (5$ is continuous at the surface of the star 

C = Co. 

Following the analysis of the gravitational potential 
equation, we introduce a second system of coordinates 
(C, /i, tp) which allo w the equation for the hydrodynamic 
potential, Eq. ( |3.6| ), to be written in a convenient form. 
These coordinates are related but not identical to those 
used by Bryan ||: 



= by/(l -e)(l-V 2 ) COS 



(4.9) 



y 



i 



(4.5) 



Thus the constant fi surfaces are hyperbolas. The coordi- 
nate fj, is confined to the range — 1 < /j < 1, with fi 2 = 1 
corresponding to the rotation axis of the star, and fi = 
to the portion of the equatorial plane outside the disk of 
radius a. The coordinate <pj measures angles about the 
rotation axis. 



The Equation (3.7) for the gravitational potential in 
these spheroidal coordinates becomes 



0_ 



(C 2 + i) 



d5$ 



d_ 



+ 



(1-M 2 ) 

C 2 + /i 2 



d 2 6$ 



(C 2 + l)(l-/i 2 ) dip 2 



0, (4.6) 



except on the surface of the star C = Co- The sepa- 
rable solutions to these equations are functions of the 
form P / m «)P, m (/x)e im¥ ' and QJ"(iC)P / m (/i)e im¥ '. The as- 
sociated Legendre functions Q™(/i) diverge at fi 2 = 1, 
consequently only the functions P/ m (/x) appear in these 
solutions. The associated Legendre functions P™(iC) di- 
verge as C' while the Q^iiQ vanish as C~^ +1 - ) in the limit 
C — > oo. Thus the gravitational potential in the exterior 
of the star, C ^ Co, must have the form: 



(5$ 



wo 



(4.7) 



for some constant a. In the interior of the star C < 
Co the situation is more complicated. Both P ; m (iC) and 



Sin <y5, 



z — o^/i . 

K 

Here the focal radius b is related to a by 

,2 



a 



4-K 2 



mi+o-K 2 



(4.10) 



(4.11) 



(4.12) 



The parameter b is real for frequencies, such as those of 
the r— modes, which satisfy k 2 < 4. As we shall see, the 



Eq. ( |3.6| ) for the potential 5U separates nicely in terms 
of these coordinates. However these new coordinates are 
rather unusual, so we present an in depth discussion of 
them in the Appendix. In summary, the coordinates £ 
and fl cover the interior of the star when their values are 
confined to the domains £ < £- < 1 and —Co < A < Co, 
where Co is defined as 



« 2 C 2 



£ 2 



2^2 



C 2 « 



b 2 4 



4(1 + Co 2 ) ^ * 



< — < 1. (4.13) 



The surface £- = 1 corresponds to the rotation axis of the 
star, and the surface fi = corresponds to the equatorial 
plane of the star. The surface of the star, C = Co, is di- 
vided into three regions in this coordinate system. The 
portions of the stellar surface nearest the two branches 
of the rotation axis correspond to the surfaces jl = ±£ D , 
while the portion of the stellar surface that includes the 
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equator corresponds to the surface £ = £ G . The coordi- 
nate p coincides with the value of the coordinate fi in 
that portion of the surface of the star where £ = £ Q . In 
the other portions of the surface of the star the value 
of fi coincides with ±£. These facts will be essential in 
imposing the bou nda ry conditions in the next section. 

The Equation (3.6) for the potential SU reduces to the 
following in terms of the coordinates (£, p, ip) 



d_ 



(e i) 



dSU 



3£ 



d_ 

dp 



(1 



p 2 ) 



p 2 



dSU 



dp 



d 2 5U 



(£2_1)(1_ M 2) d(p 2 



= 0. (4.14) 



The separated solutions of this equation are associated 
Legendre functions of £ and p. The coordinate £ in- 
cludes £ = 1 in its r ange, so the non-singular sepa- 
rated solutions to Eq. fl4.14| ) are P/ 7l (£)P J m (/i)e im¥ ' and 



Pl n (0Qr(P') eimv - The "angular" coordinate p does not 
include ±1 in its range. Thus at this stage, it is not pos- 
sible to eliminate the Ql n (p) solution without imposing 
the boundary conditions. 



V. IMPOSING THE BOUNDARY CONDITIONS 

In order to obtain the physical solutions to the stellar 
pulsation equatio ns, w e mus t now impose the boundary 
conditions, Eqs. ( |3. 10 ) and ( 3.11 ). The simplest bound- 
ary cond ition is the one that involves the derivatives of 
5$, Eq. (3.11). In the coor dina tes (fi , £, y >) used to find 
the solution for 5$ in Eqs. (4/7) and (|4.8| ), the unit nor- 
mal vector to the surface of the spheroid n a has only one 



nonvanishing component, 



i Co 2 + i 



C 2 +M 2 



(5.1) 



Thus, the normal derivatives that appear in the boundary 
condition, n a V <5$, can be expressed simply as rp'd^SQ. 
The gradient of the pressure, n a V a p that appears in 
Eq. (3.11) is given by Eq. (2.10). When evaluated at 



the surface of the spheroid this reduces to 
n a V a p = -4^GAW(C 2 + 1)(C 2 +M 2 



x(l-C cot- 1 Co). (5.2) 



Thus, using Eqs. Q . (gg), ( p] ) and Q the bound- 
ary condition Eq. (3.11) on <5$ is equivalent to 



Co + 1 <jQP(jCo) pm./ \ imip 

C 2 + l dP{"(iQ 



a 



priKo) dc 

aPT(fi)e im ' f ' + 8U 
Co(l- Co cot" 1 Co) 



P^{fi)e lm ^ 



(5.3) 



The first immediate consequence of this boundary con- 
dition is that the potential SU must be proportional to 
P; m (/Lt) on the surface of the star. In the last section we 
found that the potential SU was some linear combina- 
tion of PJ n {^)P{ n (p)e im V and P, m (£)Qf \p)e mip . As we 
show in the Appendix, the surface of the star is some- 
what complicated in the (£, p, ip) coordinate system. For 
the portion of the surface of the star that includes the 
equator, we found that £ = £ Q and p = fi. This fixes the 
angular dependence of SU. Therefore, throughout the 
star SU must have the form 



SU = (3 



P[ n (0 



P{ n {p)e im ' f 



(5.4) 



where j3 is an arbitrary constant. On the portion of the 
surface of the star that includes the equator, this ex- 
pression reduces to SU = f3P l m (p)e lmip = /?P J m (/i)e im *'. 
On the portion of the surface that includes the positive 
rotation axis, p = £ Q , the expression for SU reduces to 
SU = /3P, m (£)e""^ = (3Pl n {fi)e vmip since £ = // here. 
Finally, on the portion of the surface that includes the 
negative rotation axis, p = — £ OJ the expression for SU 
also reduces to SU = /3P ; m (-£)e m ^ = fiP[ n {fi)e lmip since 
£ = — fi here. Consequently, the potential SU reduces 
to the expression SU = f3P[ a (fi)e lm '' p everywhere on the 
surface of the star. Thus, the boundary condition on (5<I> 
reduces to 



c 2 + i rfQr«o) 

QTiKo) d( 



c 2 + i dpr(i( ) 



Co(l-CoCOt- 1 Co)' 



(5.5) 



We now see why it was necessary to obtain the so- 
lutions for SU in the strange and complicated (£,p,tp) 
coordinate system. These unusual coordinates have two 
essential properties: first they allow the solutions to 
Eq. ( [3.6| ) to be found in separated form, and second they 
have the property that one of the coordinates reduces 
on the surface of the star to the angular coordinate fi. 
This last property was needed to allow us to satisfy the 
boundary conditions using simple separated solutions for 
both SU and <5<E>. 



The boundary condition, Eq. ( 3.10 ), on the potential 
SU is unfortunately somewhat more complicated. The 
tensor Q ab that is used in Eq. ( 3.1C ) is most simply ex- 
pressed in cylindrical coordinates (see Eq. 3.3). Therefore 
it is simplest to consider the boundary conditions on SU 



in these coordinates. Let w 2 



y 2 denote the cylin- 



dric al rad ial coordinate. Then, the boundary condition 
Eq. ( 3. 10| ) can be expressed as 



(k 2 - A)n z d z 5U+K 2 n^d^5U + *^n™SU 

w 

2/ 2 6U + S<I> 

-k 2 (k 2 - 4)frp — — — . 



(5.6) 



The components of the unit normal vector to the sur- 
face of the spheroid, n a , that appear in Eq. (|5.6D can 
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be obtained by taking the gradient of the function that 
appears on the left side of Eq. (2.4): 



/ C 2 (l-M 2 ) 

C 2 + m 2 : 



(5.7) 



where A(k, ( a ) and B(( a ) are defined by 

C 2 + l dPp^o) 2m( 



A(kXo) 



^ m (UV4(C 2 + l)-« 2 d Z 



4-k 2 



(5.14) 



1 Q + 1 
C 2 + m 2 



(5. 



The part ial derivatives d z 5U and d^SU that appear in 
Eq. ( |5.6|) are more difficult to evaluate. To do this we 
must evaluate the Jacobian matrix that determines the 



coord inate transformation defined in Eqs. (4J3) through 
(4.11). The neede d par tial derivat ives are given in the 
Appendix as Eqs. ( A13 ) through ( |A16| ). These expres- 
sions can now be used to transform the derivatives d z and 
d^j needed in the boundary condition into expressions in 
terms of and d^. When evaluated on the portion of 
the surface of the star w ith £ = £ c an d jl = [i (using 



the relationships in Eqs. 4.12 and 4.13 ), we obtain the 
following 



(k 2 - 4)n z d z 6U + K 2 n m d w 8U 
k(4-k 2 ) 



C 2 + l 



<V4(C 2 + 1) - V Co 2 + ^ 



d s SU. (5.9) 



Similarly, when evaluated on the portions of the surface 
with jl — ±£ and £- = ±/x we find 



{k 2 - A)n z d z SU + K 2 n m d^SU 
k(4- k 2 ) 



V4(C 2 + 1) - « 2 



C 2 + m 2 ^ 



dj.SU. (5.10) 



Fina lly the n, we may combine the results of Eqs. (4.7) 
( |5.4| ), ( |5.9| ), and (5.10) to obtain the following represen- 
tation of the boundary condition Eq. (5.6): 



C 2 + i 



dPT(So) a 2mCo 



2(1 -Co cot" 1 Co) 



4- 



(5.11) 



We point out that Eq. ( 5.11 ) is valid on the entire sur- 
face of the star. For the portion of the surface where 
fx = — £ and £- = — ll on the surface, it is helpful to re- 
member that P™(-£ Q )P™(-pi) = F™(C )P™(^), while 

p^(-m) dPr(-to)/dt = -pr(v) dpr&ydz. 



In summ ary then the boundary conditions, Eqs. (5.5) 
and (5.11) are given by 



Co(l-CoCOt- 1 Co)' 



(5.12) 



A{nXo)8- 



2(1 - Co cot" 1 Co) 



«(a + /3), (5.13) 



B(Co 



c 2 + i dQr(iCo) Co 2 + i dpr(i( ) 



l (i(o) d( 



pr(Ko) dc 



(5.15) 



Note that £ G that appears on the right side of Eq. (|5.14| ) 
is an implicit function of k and Co as given in Eq. (|4. 13[) . 



VI. THE EIGENVALUES 



The boundary conditions Eqs. ( [5.12 ) and (5.13) can 
be satisfied for only certain values of the eigenvalue n. It 
is easy to see that the necessary and sufficient condition 
that there exists a solution to the boundary conditions is 

Q = F(k,(o) = A( K ,(o)B(Co) + 

2A(k, Co) - <oB(Co) [(1 + 3C 2 ) cot" 1 Co - 3Co] 



2Co(l-CoCOt- 1 Co) 



(6.1) 



We have verified that Eq. (|6.l|) is exactly equivalent (i.e. 
up to changes in notation) to Eq. (60) of Bryan 

We now wish to evaluate this eigenvalue equation an- 
alytically for the generalized r— modes of slowly rotating 
stars. The parameter Co determines the angular velocity 
of the star through Eq. fl2.3p . Large values of Co corre- 
spond to small angular velocities, and so we may expand 
our equations in inverse powers of Co- The leading order 
terms in such an expression for the angular velocity are 



n 2 



nGp 15C 2 



4 + o(C 



(6.2) 



We now wish to obtain solutions to the eigenvalue 
equation F(k, Co) = for large values of Co- We do this by 
expanding the expressions on the right side of Eq. (6.1) 
in inverse powers of Co^ 



F(«,c ) = -a-i)Co 



i 



dPT(K/2) 



pr(K/2) ^ 

x[l + 0(C o - 2 )]. 



Am 

\~K 2 

(6.3) 



Setting the coefficient of this lowest-order term to zero, 
we obtain an equation for the lowest-order expression for 
the eigenvalue, k : 



dP f m ( Ko /2) 
d£ 



4m 



P/>°/2). 



(6.4) 







TABLE I. The eigenvalues re D of the r-modes of the 
Maclaurin spheroids in the limit of low angular velocities. The 
frequencies of these modes are related to k by u) = (n — m)f2 
in the low angular velocity limit. Those frequencies denoted 
with * have the property that u>(ui + mtt) < 0. These * modes 
are subject to a gravitational radiation driven secular insta- 
bility. 





m = 


m — 1 


m = 2 


m = 3 


m = 4 


1 = 2 


0.0000 


1.0000 








I = 3 


0.8944 


1.5099 


0.6667* 








-0.8944 


-0.1766 








I = 4 


1.3093 


1.7080 


1.2319* 


0.5000* 






0.0000 


0.6120* 


-0.2319 








-1.3093 


-0.8200 








I = 5 


1.5301 


1.8060 


1.4964* 


1.0532* 


0.4000* 




0.5705 


1.0456 


0.4669* 


-0.2532 






-0.5795 


-0.0682 


-0.7633 








-1.5301 


-1.1834 








I = 6 


1.6604 


1.8617 


1.6434* 


1.3402* 


0.9279* 




0.9377 


1.3061 


0.8842* 


0.3779* 


-0.2613 




0.0000 


0.4405* 


-0.1018 


-0.7181 






-0.9377 


-0.5373 


-1.0926 








-1.6604 


-1.4042 








Using 


the Rodrii 


mes formula for the associated Legendre 



functions, this equation can be transformed (for the m > 
modes) into the form 



V 2 



1 



d m+1 Pi{n /2) 



ro+l 



(6.5) 



wher e P; is the Legendre polynomial of degree I. Equa- 
tion (6.5) is equivalent to Bryan's Eq. (83) H. This equa- 
tion admits I — m (or I — 1 for the m = case) distinct 
roots all of which lie in the interval —2 < n < 2 [13[ . 
For the case I = m + 1 the single root of Eq. ( |6J5 ) is 
k = 2/(m + 1). This agrees with the frequency of the 
classical r-mode of order m as found for example by Pa- 
paloizou and Pringle [||. The modes with m < are 
equi valen t to those with m > 0: if solution to 

Eq. (6.4) for some m, then — k is a solution for — m . 
In Table [j] we present numerical solutions to Eq. (6.4) 
for a range of different values of I and m. We see that 
for each value of m there exist solutions of this equation 
for each value of I > m + 1. For each value of I and m 
there are I — m different solutions. Thus, there exist a 
vast number of modes whose frequencies vanish linearly 
with the angular velocity of the star. We indicate with a 
* those frequencies in Table | that satisfy the condition 
lu(uj + mSY) < 0. The modes satisfying this condition 
would be driven unstable by gravitational radiation re- 
action in the absence of internal fluid dissipation (i.e. 
viscosity) ||. 

Next we wish to extend the formula for the eigenvalue 
to higher angular velocity. We define the next term in 
the expansion of the eigenvalue as 



+ K2 (- 2 + o((- 4 ). 



(6.6) 



Using this definition and Eq. ( |6.4| ) for n we find the next 
order term in F(k, £„) to be 



F(K, Co) = -(/-!) 



m 

T 



Ko l(l + l) 2k 2 1{1 + 1) 



+ §(2Z + l) Ko + 0(C 2 ). 



(6.7) 



We can now determine the second order correction to the 
eigenvalue of the mode by solving F(k, Co) = for k 2 : 



k 2 = 



«o(4 



21(1+1) 



rn 
2k 



1(1 + 1) 2 21 + 1 



5 1-1 



In the case of the classical r-modes with I = m 
expression reduces to: 



4m 



n 2 



iy 



(m+l) 2 



2 2m + 3 
5 m 



(6.8) 
1, this 

(6.9) 



Thus, for the classical r— modes the frequency of the 
mode is given by |14 



u> + mfi = fi< 



1 

3m 



(m + 1) 

l5\ 



2m + 3 



m(m + l) 5 



-0(f7 5 



n 2 I 

irGpf 
(6.10) 



Fo r ra pidly rotating stars we must solve the eigenvalue 
Eq. (|6.l|) numerically. This presents certain unusual nu- 
merical difficulties. In particular the associated Leg endre 
functions P/"(zC) and Q^^O that appear in Eq. ( 5.15| ) 
are not commonly encountered. Thus, we digress briefly 
here to describe how these functions may be evaluated 
numerically. First we note that these functions are essen- 
tially real. In particular the functions P™(C) and <5™(C) 
defined by 



pr(i() = i l Pr(o, 



0T(C), 



(6.11) 



3.12) 



are real for real values of (. The functions P™(C) can 
be evaluated numerically using essentially the same al- 
gorithms used to evaluate their counterparts on the real 
axis. In particular 



P™(C) = (2m-l)!!(C 2 + ir/ 2 , 



Pr +1 (C) = (2m + l)CP™(C). 



(6.13) 



(6.14) 



These expressions can be used as initial values for the 
recursion formula, 

(i - m)Pr(o - (2i - i)cp ; m i(o + (i+m- i)p ; m 2 (o, 

(6.15) 
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FIG. 1. Angular velocity dependence of the eigenvalues k 
of the classical r— modes, i.e. those modes with I = m + 1 for 
/ < 6.. The frequencies of these modes are related to k by 
uj — (k — m)0,. 
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FIG. 2. Angular velocity dependence of the frequency 
= (k — m)f2 of the classical m = 2 r— mode. The solid curve 



gives w corresponding to the exact solution of Eq. (3.1), while 
the dot-dashed and dashed curves correspond to the first and 
second order approximations respectively from Eq. (6.10). 



from which the higher order functions can be determined. 
This approach does not work for the Qf 1 however. The 
problem is that we need to evaluate these functions over a 
fairly wide range of their arguments, e.g. 0.25 < ( < 75. 
For large values of ( the recursion for the involves 
a high degree of cancellation among the various terms. 
Standard computers simply do not have the numerical 
precision to perform these calculations to sufficient ac- 
curacy. Instead we rely on an integral representation of 
Qf - Based on Bateman's Eq. (3.7.5) @, we find 



l (C) 



(-1) 



l+m+1 



(I 



m 



2'+ 1 Z!(C 2 + l) m / 2 

a2 _|_ £2^(rn-;-l)/2 

-1 

X COS 



X 



(l-t 2 ) 1 



[(I + 1 -m) tan' 1 (t/£)]dt. (6.16) 

The integrand in this expression is well behaved for all 
values of £ > 0, and the integrals may be determined 
numerically quite easily. 

Using these numerical techniques, then, it is straight- 
forward to solve for the eigenvalues of Eq. (3.1), 
F(k, £ ) = 0, over the relevant range of the parameter 
0.25 < Co < 75. Figures [l], |], and || display the angular 
velocity dependence of the eigenvalue k for a number of 
r— modes. Figure |l| depicts k for the 'classical' r— modes, 
I — r7i+l, with I < 6. This / = 777+ 1 = 3 mode is the one 
found by Lindblom, Owen, and Morsink || to be suffi- 
ciently unstable due to gravitational radiation emission 
that it is expected to cause all hot young rapidly rotat- 
ing neutron stars to spin down to low angular velocities 
within about one yea r. W e have verified that our numer- 
ical solutions of Eq. (|6.l|) agree with those of Eq. ( 3.10 ), 
up to terms that scale as f2 5 . 

Figure |^ presents the angular velocity dependence of 
the frequency of the classical I — m + 1 = 3 r— mode 
as measured in an inertial frame. The solid curve corre- 
sponds to the exact solution to Eq. (|6.1|), while the dot- 
dashed and dashed curves represent the first and second 



order approximations (respectively) given in Eq. ( 3.10 ). 
The units in F ig. [2| for both the vertical and horizontal 
axis scale as ■\pKGp. The value of p = 7 x 10 14 gm/cm 3 
chosen here represents a typical average density for a 
neutron star. Figure || illustrates three interesting fea- 
tures about the frequencies of this mode in large angu- 
lar velocity stars. First, the frequency is only about 2/3 
that predicted by the first-order formula for stars rotating 
near their maximum angular velocity. This means that 
the gravitational radiation reaction force, which scales as 
(u} + mQ)uj 5 , could be about 1/5 of that predicted by the 
first order formula. (Unless the mass and current mul- 
tipoles of the rapidly rotating models are much larger 
than their non- rotating values.) Second, the frequencies 
of these modes first increase and then decrease as the an- 
gular velocity of the star is reduced. This means that the 
time evolution of the gravitational radiation signal from 
these sources will be more complicated than had been 
anticipated by Owen, et al. M on the basis of the first or- 
der expression for the frequency. In particular it appears 
that the evolution of the frequency will not be monotonic 
for the most rapidly rotating stars. Third, the accuracy 
of the second-order formula for the frequency is in fact 
considerably worse than that of the simple first-order for- 
mula for very rapidly rotating stars. This suggests that 
any application of low angular velocity expansions of the 
r— modes to the study of rapidly rotating stars is some- 
what suspect. 

Figure || depicts the angular velocity dependence of k 
for several other previously unstudied r— modes. The 
modes depicted in Fig. || all have the property that 
lu(oj + mO) < 0, and hence these modes would all be 
subject to the gravitational radiation secular instability 
in the absence of internal fluid dissipation (i.e. viscos- 
ity). These additional modes all couple to higher order 
gravitational moments than the classical I = m + 1 = 3 
mode. Thus, these additional modes probably do not 
play a significant role in the astrophysical process which 
spins down hot young neutron stars. 
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FIG. 3. Angular velocity dependence of the eigenvalues k 
of r— modes with 5 > I > m + 1 which are unstable to the 
gravitational radiation instability. The frequencies of these 
modes are related to n by uj — (k — m)Q. 

VII. THE EIGENFUNCTIONS 

The eigenfunctions associated with these modes are 
those given in Eqs. (4.7), (4.S), and (5.4). Inside the 
fluid we have 



i \ l Qo) 



~\ imtp i(K-m,)Qt 



while outside we have 



5$ 



a 



WO 



iT(M)e 



(7.1) 



(7.2) 



(7.3) 



The coefficients a and (3 are determined by solv ing the 
boundary conditions. The eigenvalue Eq. (6.1) is the 
consistency condition for the existence of these solutions. 
Given a solution of the eigenvalue equatio n the n, the ratio 
of a and /3 can be determined from Eq. ( 5.1 2| ) : 



P = - [1 + Co(l - Co Cot" 1 Co)-B(Co)] O- 



(7.4) 



It might appear at first glance that the eigenfunctions 
associated with the I — m distinct r— modes are identi- 
cal. However, the coordinates £ and fx depend on the 
eigenvalue k. Thus, the spatial dependence of SU will be 
different for each of these modes. Interestingly enough, 
however, the spatial dependence of the gravitational po- 
tential, 5$, depends only on / and m, and hence is the 
same for alH — m distinct modes. 

While the expressions for the eigenfunctions are quite 
simple in terms of the special spheroidal coordinates used 
here, they are rather complicated when expressed in more 
traditional coordinates. One exception to this is the case 
of the classical r— modes, i.e. those with / = m + 1. In 
this case the needed associated Legendre function, -P™ +1 , 
has the simple expression 



^Vi(C) = (-l) m (2m + 1)!!C(1 - C 2 ) m/2 . (7.5) 

Using the expressions for the bi-spheroidal coordinates 
(£, /Et) in term s of t he c ylindrical coordinates (w, z) given 
in Eqs. and (|A7|), it follows that 



^i(O^h-iCA) 



(-ir(2m+l)!!P- +1 (£ ) 



Thus, the hydrodynamic potential SU is given by 



SU = f3—— ( — ) e imtp e i(K-m)nt_ 



(7.6) 



(7.7) 



Aside from the overall normalization then, the spatial de- 
pendence of SU is completely independent of the angular 
velocity of the star! A similar expression can be obtained 
for the gravitational potential (5$ within the star. Using 
the fact that 

= (-l) m (2m + l)!!iC + i(iCo)-f f- , (7.8) 

tip V-fle/ 



the gravitational potential is given by 



5$- 



Rn 



e imip gi(K- 



(7.9) 



Thus the spatial dependencies of the potentials SU and 
SQ are identical for the classical r— modes! This spatial 
dependence can also be expressed in spherical coordinates 
as r m+1 Y m+ i m (9, </?), up to an overall normalization. It 
is interesting that this same function satisfies both the 
hydr ody namic Eq. (|]^) and the gravitational potential 
Eq. (@. 

The Eulerian perturbation in the pressure for these 
modes is determined from these two potentials by 
Eq. (O: 



(7.10) 



= (^ a + f3)— ( — ] ,.'»',-, 'I'.' ""•-'-i 



R v V R 



We note that the constants a and j3 used here have been 
scaled by the factor (— l) m (2m + 1)!! com pared to their 
original definitions in Eqs. (4.8) and (5.4). It is also in- 
structive to evaluate the Lagrangian perturbation of the 
pressure, Ap, as d efine d in Eq. (pjj). Using the expres- 
sions in Eqs. fl7.7| ) and (|7.9D, we find that 




2(1 - Co cof 1 Co)[(2 - «)(! + C 2 ) ~ ™<l 
K 2 (2 - k)C„[(1 + 3C 2 ) cot" 1 Co - 3Co] 



>SU 



(7.11) 
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The term enclosed in {} brackets in Eq. (7.11) depends 
only on the frequency of the mode k, the angular ve- 
locity of the star (through ( a ), and the amplitudes of 
the perturb ation s a and /3. Wh en th e boundary con- 
dition Eq. (3.10) (or equivalently 5.13) is satisfied, this 
term vanishes. Thus, we find that the Lagrangian per- 
turbation in the pressure Ap vanishes identically for the 
classical r— modes of the Maclaurin spheroids. And this 
result (which is a consequence of the extremely simple 
eigenfunctions for these modes) holds for stars with any 
angular velocity! 



VIII. DISCUSSION 

Our analysis, which follows closely in the footsteps of 
the remarkable analysis of Bryan, provides several in- 
teresting insights into the properties of the r— modes 
of rapidly rotating stars. It demonstrates for example, 
that these modes actually do exist in rapidly rotating 
barotropic stars, and are not just figments of the first- 
order perturbation theory (as claimed by some authors 
p5|). This analysis shows that some properties of the 
r— modes are not well approximated by the low angular 
velocity expansions. Figure [2] illustrates, for example, 
that the first-order expression for the angular velocity 
dependence of the frequency of the classical r— modes 
is in fact superior to the second-order expression in the 
most rapidly rotating stars. This analysis shows that the 
frequency evolution of the gravitational radiation emit- 
ted by the r— mode instability is likely to be more inter- 
esting than had previously been thought Qj. Figure || 
shows that the frequency of these modes will first in- 
crease and then decrease as the angular velocity of the 
star is reduced by the emission of gravitational radia- 
tion. This analysis shows that there is a much larger 
family of r— modes than the 'classical' r— modes studied 
for example by Papaloizou and Pringle Q . For each pair 
of integers I and m (which satisfy I > m > 0) there exist 
I — m (or I — 1 in the m = case) distinct r— modes. 
This analysis shows that a significant fraction of these 
previously unstudied r— modes are subject to the gravita- 
tional radiation driven secular instability. This analysis 
has derived simple analytical expressions for the eigen- 
functions of the classical r— modes. Both of the poten- 
tials SU and <5$ are proportional to r m+1 Y m+lm (0, <p) for 
the classical r— modes in Maclaurin spheroids of any an- 
gular velocity. This analysis shows that the Lagrangian 
variation in the pressure, Ap, associated with the classi- 
cal r— modes vanishes identically in Maclaurin spheroids 
of arbitrary angular velocity. Thus, the r— modes of 
the Maclaurin spheroids provide a completely unsuitable 
model for the study of the effects of bulk viscosity on 
the r— modes. This analysis also provides an interesting 
mathematical example of a hyperbolic eigenvalue prob- 
lem. The r— modes studied here have the property that 
k 2 < 4. Thus, the equation satisfied by the potential SU, 



Eq. ( p.6[ ), is in fact hyperbolic. Nevertheless, the b ound - 
ary condition imposed on the potential SU, Eq. ( 3.10 ) 



is of the mixed Dirchlet-Neumann type that is generally 
associated with elliptic problems. The analytical solu- 
tions given here illustrate that this unusual hyperbolic 
eigenvalue problem does nevertheless admit well behaved 
solutions. 
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APPENDIX: BI-SPHEROIDAL COORDINATES 



Th e coordinates £ and jl defined in Eqs. (4.9) through 
(4.11) are rather unusual. The purpose of this appendix 
is to explore the geometrical properties of these coordi- 
nates. The coordinates £ and jl cover the planes usually 
described with the spherical coordinates r and 9 or equiv- 
alently the cylindrical coordinates va — y x 2 + y 2 and z. 
First, we note that it follows from Eqs. (|4.9|) through 



(4.11) that 



— - 



(Al) 



Thus, the surfaces of constant £ (with £ 2 < 1) are 
spheroids for the r— modes which have k 2 < 4. The par- 
ticular surface £ — £ , with £ Q defined by 



? 2 



a 2 Q 



a* 2 



b 2 4-K 2 4(l + Cf)-'« 1 



< T <1 > 



(A2) 



is identical to the surface of the star £ = ( . Paradoxi- 
cally, the constant jl surfaces also satisfy the equation 



—- 



1 - ji 2 



4 — k 2 jl 2 



(A3) 



Thus, the constant jl surfaces are the same family of 
spheroids as the constant £ surfaces. Thus the coordi- 
nates £ and jl constitute a bi-spheroidal coordinate sys- 
tem. We note that the equatorial and polar ra dii o f the 
spheroid, R e and R p defined in Eqs. ( gjj ) and (2.6), are 
related to the constants b and £ by 



R 2 =b 2 (i-e ), 



(A4) 
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(A5) 



In order to understand how the coordinates £ and p 
cover the interior of the star it is helpful to introduce 
two additional coordinates s and 9: 



s 2 sin 2 I 



w 



i-C 2 



s 2 cos 2 1 



2 7,2 



R 2 



(A6) 



(A7) 



The coordinate s has the value 1 on the surface of the 
star and at its center. The coordinate 9 ranges from 
the value on the positive rotation axis, through ir/2 
on the equatorial plane, to it on the negative rotation 
axis. Thus, the coordinates s and 9 map the interior of 
the star into the unit sphere in a natural way. It will 
be instructive to express the coordinates £ and p then 
in terms of s and 9. These expressions can be obtained 
from Eqs. (TaI) and (|A7|): 



( 2 = k(u + v), 



(A8) 



(1 + s(, cos + s-y/l - H si n 9) 
x (1 + s£ cos9 — syl — £ 2 sin( 
x (1 ~ s £o cos + syl — £ 2 sin ( 
x (1 + s£ D cos# — syl — £ 2 sin( 



(A12) 



Each of the terms on the right side of Eq. ( |A12| ) is pos- 
itive, since each is 1 plus the inner product of a pair of 
unit vectors multiplied by s. Thus v 2 > and so v is 
real and positive. Further, v < u and so £ 2 and p 2 are 
positive. Thus £ and p are finite and real for each point 
of the interior of the star. 

Second we wish to show that the transformation be- 
tween the bi-spheroidal coordinates (£, fl) and the stan- 
dard cylindrical coordinates (m, z) is non-singular every- 
where within the star. We do this by evaluating the 
Jacobian matrix (i.e. the matrix of partial derivatives) 
of the transformation: 



dz 



Kt(i-e) 



b(jP-?)VF- 



dc _ - em - a 2 ) 

dm b(fl 2 - £ 2 ) 



(A13) 



(A14) 



where 



fl 2 = \{u- 



(A9) 



dp 

dz 



<(i - p 2 ) 



b(p 2 -£ 2 )V4^? 



(A15) 



l-s 2 + S 2 (C 2 + > 



4s 2 £ 2 c °s^ 0. 



(A10) 



(All) 



We note that while Eqs. (A6) and (A7) are symmetric 
in £ and fl, this symm etry has b een broken in order to 
obtain the expressions (AS) and flA9). 

We will now show that the coordinates £ and fl also 
cover the interior of the star when their values are re- 
stricted to the ranges: £ Q < £ < 1 and — £ D < fl < £ D . It 
is easy to see from the second equalities in Eqs. ( |A6| ) and 
( A7) that each point (£, /i) in the domain £ Q < £ < 1 and 
— Co < M < Co corresponds to a point within the star (i.e. 
a point with s < 1). Proving the converse is more diffi- 
cult. We do this in three steps: First we show that the 
functions £ and fl are real and finite at each point within 
the interior of the star. Second we show that these func- 
tions have no critical points (e.g. no maxima or minima) 
except on the surface of the star. Third, and last, we 
show that £ and fl are confined to the ranges £ Q < £ < 1 
and — £ D < ft < £ for points on the boundary. 

First we show that £ and fl are real and finite for 
each point within the star. The quantity u defined in 



Eq. ( A10) is positive for points within the star (i.e. points 
with s < 1). Next we show that v (defined as the posi- 
tive root in Eq. All) is real and thus positive for points 



dp /V(i-£ 2 )(i-ft 2 ) 

dm KP 2 -e) 



(A16) 



within the star. We do this by re-writing v as 



These expressions show that the Jacobian matrix, and 
hence the coordinate transformation between (£, p) and 
(m,z), is non-singular except for the points within the 
star where £ = 1, or £ 2 = p 2 — £ 2 . The transformation is 
also singular at p 2 = 1, however, these points are not in 
the range of interest to us here. The non-singularity of 
the Jacobian matrix proves that V a £ and V a /2 are non- 
vanishing everywhere within the star. Thus, the maxi- 
mum and minimum values of £ and p will only occur on 
the surface or the rotation axis of the star. 

Third, and finally, we explore the values of the coordi- 
nates (£, p) at specific physical locations in the star, e.g. 
the surface of the star, the rotation axis, etc. We be- 
gin first with th e rot ation axis . In the (s, 9) coordinates 
defined in Eqs. (A6) and (A7), the rotation axis corre- 
sponds to the points where sm.9 = 0. If follows then that 
these points correspond to £ = 1 and p 2 = s 2 £ 2 . Thus, 
the rotation axis is the surface £ = 1, a singular surface 
of the coordinate transformation. The coordinate p takes 
on its entire range, — £ e < p < £ OJ for points along this 
axis. The equatorial plane of the star, cos 9 = 0, corre- 
sponds to the coordinate surface p = 0. The coordinate 
£ 2 = 1 — s 2 + s 2 £ 2 ranges from £ = 1 (on the rotation 
axis) to the value £ = £ Q on the surface of the star. 
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The surface of the star, s = 1, is unexpectedly compli- 
cated in the (£, p) coordinate system. For points of the 
stellar surface near the equator, cos 2 9 < £ 2 , the function 
v defined in Eq. (All) has the value: v — £ 2 — cos 2 9. 
Thus, the surface of the star in this region has £ = £ Q 
while p — cos 9. In the regions of the stellar surface 
near the rotation axis, cos 2 9 > the function v (de- 
fined as the positive root in Eq. All) has the value: 
v — cos 2 9 — £ 2 . Thus the stellar surface in these regions 
have p = ±£ Q and £ 2 = cos 2 9. The role of £ and p as 
"radial" and "angular" coordinates are reversed therefore 
in different regions of the star. 

In summary then, we have shown that the maximum 
and minimum values of £ are 1 and £ respectively, while 
the maximum and minimum values of p are ±£ D . This 
concludes our demonstration that the coordinates (£, p) 
when confined to the ranges £ < £ < 1 and — £ c < p < (,o 
do form a non-singular coordinate system that covers the 
interior of the star. The transformation between these 
and the usual cylindrical coordinates is singular only on 
the rotation axis, £ = 1, and at the singular point £ 2 = 
P- 2 = £o on the surface of the star. 

Finally, it will be useful to work out the relationship 
between the surface values of the (£, p) coordinates with 
those of the oblate spheroidal coordinates (C, A 4 )- Us- 
ing the defi nition s of the (s,9) coordinates introduced 
in Eqs. (A6) and (A7), and the definitio ns o f t he o blate 
s phe roidal coordinates (£, //) from Eqs. (4.1), (4.2), and 
(4.3), it is straightforward to show that 



s 2 sin 2 9 



(C 2 + i)(i-M 2 ) 

Q + 1 



s 2 cos 2 i 



cv 

C 2 



Thus on the surface of the star, £ = C and s — 1, we 
have cos# = /i. We have also shown that p = cos 9 on 
the portion of the surface of the star where cos 2 9 < £ 2 , 
and that £ = ± cos 9 on the portion of the surface where 
cos 2 9 > £ 2 . Thus we find that p = \i on the portion of 
the surface where £ = £ , and that £ = ±/z on the portion 
of the surface where p — ±£ Q . 



N. Andersson, Astrophys. J. , in press, gr-qc/9706075 
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